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Abstract 

This study is concerned with how the attractor dimension of the two-dimensional 
Navier-Stokes equations depends on characteristic length scales, including the sys- 
tem integral length scale, the forcing length scale, and the dissipation length scale. 
Upper bounds on the attractor dimension derived by Constantin-Foias-Temam are 
analysed. It is shown that the optimal attractor-dimension estimate grows linearly 
with the domain area (suggestive of extensive chaos), for a sufficiently large domain, 
if the kinematic viscosity and the amplitude and length scale of the forcing are held 
fixed. For sufficiently small domain slightly "super-extensive" estimate be- 

comes optimal. In the extensive regime, the attractor-dimension estimate is given 
by the ratio of the domain area to the square of the dissipation length scale defined, 
on physical grounds, in terms of the average rate of shear. This dissipation length 
scale (which is not necessarily the scale at which the energy or enstrophy dissipation 
takes place) can be identified with the dimension correlation length scale, the square 
of which is interpreted, according to the concept of extensive chaos, as the area of 
a subsystem with one degree of freedom. Furthermore, these length scales can be 
identified with a "minimum length scale" of the flow, which is rigorously deduced 
from the concept of determining nodes. 
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1 Introduction 



It has become a classical idea that the attractor dimension — a measure of dy- 
namical complexity — of a certain class of dissipative physical systems scales 
linearly with the volume of the system, a property known as "extensive chaos" 
[8,15,16]. This central idea of extensive chaos is physically interpreted and best 
appreciated with the following heuristic argument. In a dissipative system 
there exists a dissipation length scale £, below which all modes are damped; 
only modes with length scales above the dissipation scale are dynamically ac- 
tive. Intuitively, the attractor dimension is approximately equal to the num- 
ber of active degrees of freedom. If other parameters can somehow be con- 
trolled so that £ remains fixed while the volume V of the system of typi- 
cal length scale L is increased, then the number of active modes grows as 
(L/£) d where d is the spatial dimension of the system — i.e. the degrees of 
freedom grow linearly with V . For demonstrative purposes let us consider a 
dissipative dynamical system defined in a (Hilbert) space H(Q) of doubly- 
periodic functions on Q = (0,2irL) x (0, 2irL). This space is spanned by 
{cos L^i^kiXi + ^2^2), sin L~ l {kiX\ + ^2^2)} where ki and k 2 are integers. A 
basic mode of wavevector (k\, k 2 ) T above the length scale I satisfies 

(jfe2 + £2)1/2 < L/L 

Modes satisfying the above inequality can be geometrically identified as lattice 
sites inside the circle of radius L/£, the number of which is approximately the 
area of the circle (for a sufficiently large radius). It follows that the number 
of modes above the length scale £ is proportional to (L/£) 2 which is propor- 
tional to the domain area. Hence, extensive behaviour might be expected to 
be common among dissipative dynamical systems for which the dissipation 
length scale is independent of the system size. (Thus, the problem of proving 
that a system is extensively chaotic might be restricted to showing that its 
dissipation length scale is size- independent.) Although expected to be com- 
mon, extensive behaviour has only been partly validated for a few cases. On 
the theoretical side, for example, extensive behaviour has been suggested for 
the 2D Navier-Stokes equations on an elongated domain for a special class of 
the driving force [38] and for the one-dimensional complex Ginzburg-Landau 
equation [1]. 2 On the experimental side, it has been shown numerically that 
the Kuramoto-Sivashinsky model [24] and the Miller-Huse model [26] behave 

2 In those works, an upper bound on the attractor dimension that scales linearly 
with the domain "volume" is derived. Strictly speaking, such an upper bound does 
not fully establish extensivity, without a lower bound that behaves in the same man- 
ner. For the 2D Navier-Stokes equations on an elongated (square) doubly periodic 
domain driven by a simple forcing, a lower bound on the global attractor that scales 
linearly with the domain area is derived by Babin and Vishik [2] (Liu [23]). 



2 



extensively. 

Another physical interpretation of the idea of extensive chaos, according to 
its originators, is that large systems become complex with increasing size in 
a simple way by replicating weakly-interacting and statistically similar sub- 
systems of some characteristic size [8,15,16]. In other words, an extensively 
chaotic system is a system such that there are no interesting collective effects 
with growing size, as might be expected if small dynamical units somehow 
bind themselves into a larger effective unit with fewer total degrees of free- 
dom. The characteristic size is identified as the dimension correlation length 
[8,9,10,15,16], defined as £5 = (V/D) 1 ^, where V is the volume of the system, 
D is the (global) attractor dimension, and d is the dimension of the system's 
physical space. Hence, a subsystem with volume = V/ D is interpreted as a 
dynamical unit with one degree of freedom. Clearly, if £5 remains fixed as V is 
varied, then D oc V and the system is extensive. So in this case the problem 
of proving extensivity amounts to showing that £5 is size-independent. 

In this study we examine the extensivity of two-dimensional turbulence. The 
two standard estimates of attractor dimension derived by Const ant in-Foias- 
Temam [6,7] for the Navier-Stokes equations with doubly periodic boundary 
condition are analysed. It is shown that one of the estimates, when modified to 
take into account the forcing length scale, is optimal, from the perspective of 
extensive chaos. Moreover, the modification makes the physical significance of 
the forcing length scale more apparent. In particular, a forcing with fixed am- 
plitude, when applied at progressively smaller length scales, injects a smaller 
amount of energy, leading to an increase in the dissipation length scale (de- 
crease of Reynolds number and of attractor dimension). The other estimate, 
which is slightly super-extensive, turns out to be optimal when the system's 
length scale is sufficiently small. We also address the physical relevance of the 
dimension correlation length, and its interpretation as the size of relatively 
independent subsystems in an extensively chaotic mother system. In fact, the 
dimension correlation length can be identified with the dissipation length as 
defined, on physical grounds, by Hohenberg and Shraiman [16]. Furthermore, 
this length scale is equivalent to a "minimum length scale" rigorously deduced 
from the concept of determining nodes [13,14]. 



2 Attractor-dimension estimates for 2D turbulence 

The asymptotic behaviour of the solutions of the forced two-dimensional Navier- 
Stokes equations: 

— + (u- V)w + Vp = vAu + f(x), V-m = 0, xencTZ 2 , (1) 
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has been the subject of intense research for the past 30 years [2,5,6,7,11,12,21,33,35,36,38]. 
The existence of a bounded finite-dimensional (Hausdorff or fractal dimension) 
global attractor has been established [12,21], and estimates on its dimension 
calculated [2,6,7]. The idea behind this calculation is that if an arbitrary Tri- 
dimensional volume element of the phase space asymptotically compresses to 
zero, then the number of positive Lyapunov exponents of the system and the 
Hausdorff (as well as the fractal) dimension of the global attractor cannot 
exceed m. The greatest lower bound on all such possible m is then an up- 
per bound on the attractor dimension. 3 The book of Temam [32] provides 
details of this mathematical technique and its application to various dissi- 
pative dynamical systems. Two standard attractor-dimension estimates are 
obtained when doubly periodic solutions on Q = (0, 2nL) x (0, 2nL) of (1) are 
considered: 



D< Cl G, (2) 
D<c 2 G 2/3 (l + \nG) 1/3 , (3) 

where D is the attractor dimension, c\ and c 2 are absolute constants, and G 
is the generalized Grasshof number defined by 



G^. (4) 

In (4) ||-|| is the £ 2 -norm and Ai = L~ 2 is the smallest eigenvalue of —A in the 
solution phase space, being the space of non-divergent zero-average vector- 
valued functions on Q. These estimates rely on the asymptotic averages 



1 

\u\\i) = limsup — dt u(—A)udx, (5) 
o n 

T 

\uf 2 ) = limsup — f dt [ | Au\ 2 dx. (6) 



o n 



Note that and ^(||w|| 2 ) are, respectively, the average energy dissipation 

(hence, the average energy injection) and the average enstrophy 4 dissipation 



3 By definition, the Hausdorff dimension of a set is always bounded from above by 
its fractal dimension and the upper bound can in cases be strict. Nevertheless, the 
greatest lower bound on m derived in this manner is an upper bound on both the 
fractal and Hausdorff dimensions of the global attractor (cf. [4,27]; [31], Appendix 
B). Hence, there is no need to distinguish between the two dimensions in this context 
and the term "attractor dimension" hereafter may refer to either one. 

4 The quadratic quantity \\u\\ 2 /2 (or just ||m||i) i s known in the classical theory of 
turbulence as the enstrophy. 
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(hence, the average enstrophy injection). For fixed average forcing amplitude 
11/11 /L and fixed viscosity, G oc L 3 so (2) and (3) are not linear in the do- 
main area. Hence, these attractor-dimension estimates are not suggestive of 
extensive chaos although (3) is "essentially" extensive, only slightly "super- 
extensive" by a logarithmic term. 

However an intermediate step towards (2) is 



o<_ Cl Mp. (7) 

The energy equation gives the straightforward asymptotic inequality ( 11^ II i) ly/2 < 

1 /2 

11/11 /(^-V )) from which (7) then implies (2). As noted independently by 
Robinson [30] and Tran [34], the two attractor-dimension estimates (2) and (3) 
do not take into account the spatial scale of the forcing. A simple modification 
to the estimate of (||w||i) improves (2) and gives it an explicit dependence on 
the forcing length scale (see [2,30]; [34], p. 40). To this end let us consider the 
energy evolution equation 



2dt 



| u || 2 = -is\\u\\l + I ufdx = -u\\u\\l + J(-A) 1/2 u(-Ay 1/2 fdx 



1 d 

\\u\l" = —v \u\\ + / u j dx = —v \u\\ + 
n n 

< + INIi ll/IU<4HI< + ^II/1-i> (8) 



where the two inequalities are obtained by applying the Schwarz and Young 
inequalities, respectively. It follows that asymptotically 



^><tW|I-i- (9) 



v 



In terms of this estimate (7) becomes 



D<c^ 2 . (10) 

This estimate, together with (3), which is written explicitly in terms of the 
physical parameters as 



2/3 / II XII \ V 3 



are the focus of our present analysis. 
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The norm ||/||_ 1 implicitly contains information about the spatial scale of /. 
For example, a monoscale forcing /, i.e. —A/ = A s / where A s is an eigenvalue 
of —A, yields ||/||_ 1 = A7 1 / 2 ||/||. For a general forcing one may define the 
forcing length scale £ s = || y 1 1 x / ||/|| [34]. Note that i s is always 5 less than L, 

1 /2 

so that \\fW_! < 11/11 /A/ . Hence, (10) is an improved estimate as compared 
to (2) (see [34], p. 44). This fact is independently noted by Robinson [30]. 

We note in passing that while the energy injection decreases when a forcing 
with fixed ||/|| is applied at progressively smaller length scales, the enstrophy 
injection j n (—A)u f dx does not necessarily do so. Therefore, it does not seem 
possible to apply the above trick to improve (3), which is based on an estimate 
of the enstrophy injection. 

Remark 1. The appearance of the forcing length scale in (10) is natural. The 
forcing scale has been shown to play a central role in the spectral distribution 
of energy and enstrophy and of their dissipation [36]. It also has non-trivial 
effects on optimal conditions for nonlinear stability [37]. 

Remark 2. It is a highly non-trivial problem to determine the energy injection 
rate J n u f dx. This depends on how the fluid responds to a particular forcing / 
and cannot be identified with ||/|| 2 . The estimate of J n u / dx using the Schwarz 
and Young inequalities, leading to (9) and eventually to (10), provides what 
is essentially an upper bound on the energy injection rate. When a forcing 
/ with fixed ||/|| is applied at progressively smaller scales, this upper bound 
decreases, leading to a decrease of the attractor-dimension estimate (10). The 
interpretation is that the decrease of the attractor dimension is due to the 
decrease of the energy injection, leading to an increase in the dissipation length 
scale (decrease of Reynolds number). 

Remark 3. It is concluded in [30] that when a fixed amount of energy is 
injected at progressively smaller scales then the dimension of the attractor 
decreases. This conclusion seems to originate from the (mis) identification of 
||/|| 2 with the energy injection rate. It is correct to say that a forcing with a 
fixed [[/[I when applied at progressively smaller scales leads to a decrease of 
the attractor dimension, as seen above. However, such a forcing, when applied 
at progressively smaller scales, does not inject the same amount of energy and 
the interpretation in the preceding remark is more plausible. 

Remark 4. To clarify the arguments in the previous two remarks let us con- 
sider the case of monoscale forcing /. This forcing gives rise to the (primary) 

5 The case i s = L is pathological as the forcing then belongs to the eigenspace 
of —A corresponding to Ai. This type of forcing is well known to lead to a trivial 
attractor consisting of a single stable stationary solution. This fact was first proven 
by Iudovich [17] and later investigated by Marchioro [25]. A simple proof is given 
by Constantin-Foias-Temam [7]. 
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stationary solution u = (— z/A) -1 / = f/u\ s (see Remark 5), which is sta- 
ble for sufficiently small ||w|| [37]. In the stable regime, the energy injection 
rate is J n uf dx = ||/|| 2 /(^A S ) = ||/|| 2 _ 1 /za This energy injection rate indeed 
decreases with increasing A s (decreasing forcing length scale), for a fixed ||/||. 

Remark 5. The fact that a monoscale forcing gives rise to a monoscale 
stationary solution is more readily deduced from the vorticity equation — 
d t Aip + d x ip dyAip — dyip d x Aip = vA 2 i\) + g, where g = z ■ V x / (z is the 
normal to the fluid domain) — than from (1). If g satisfies Ag = —\ s g, then 
■0 = —g/(u\ 2 ) is a stationary solution because the nonlinear terms identi- 
cally vanish for ip — ■0. In the traditional NS system, (u ■ V)w is non-zero 
(and not even non- divergent) for a non-divergent monoscale u. However, the 
divergent component of (u ■ V)u is balanced by the pressure gradient term, 
and its non-divergent component vanishes for a non-divergent monoscale u. 
The latter component is described in the literature as the projection of the 
nonlinear term onto the non-divergent zero-average phase space, and is often 
denoted by B(u, u). This projection is essentially equivalent to taking the curl 
of (u ■ V)u to obtain the vorticity equation. Detailed treatments of B(u,u) 
are given in [23,37], where it is noted that the nonlinear interaction coeffi- 
cient between two modes k and I contains the factor z ■ k x I (\k\ 2 — |/| 2 ) [z is 
the normal to the wavevector plane]. (Pedlosky [28] notes that the nonlinear 
term in the vorticity equation vanishes for waves with the same wavelength or 
with parallel wavevectors.) Therefore, if the forcing / in (1) consists of basic 
modes, for which the wavevectors are parallel or have the same magnitude, 
then (— vA)~ l f is a stationary solution. Thus, u is a stationary solution. 



3 Characteristic length scales 

There are various length scales that characterize the dynamics of a spatially 
extended system. Hohenberg and Shraiman [16] distinguish three length scales 
which are associated with excitation, dissipation, and correlation, and suggest 
that it is the ratios of these length scales to one another and to the typical 
system length scale that will determine the dynamics of the system. The exci- 
tation length scale is straightforward; it is the scale on which energy is injected 
into the system by external forcing. The dissipation length scale I characterizes 
the length scale below which all modes are damped in a finite time. For fluid 
flow where molecular viscosity alone is responsible for dissipation, t ~ (z/r) 1 / 2 
where r -1 is the local rate of shear [16]. This definition is applicable to both 
two- and three-dimensional flows. However, it is not the only definition of I in 
each case. In two-dimensional flow, for example, Batchelor [3] and Kraichnan 
[20] define the dissipation length scale as t = (z/ 3 /^) 1 / 6 , where i] is the average 
enstrophy flux to small scales which asymptotically must equal the average 
rate of enstrophy dissipation u(\\u\\l)/L 2 . It should be noted that these defini- 
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tions have their roots entirely in physical and dimensional considerations. In 
particular, the latter definition is based on the hypothesis that the enstrophy 
cascades to small scales (< £) via a k~ 3 energy spectrum, so that £ is approxi- 
mately identified with the length scale corresponding to the high-wavenumber 
end of the enstrophy-cascading range [20] . In the absence of a direct enstrophy 
cascade on the global attractor, as is argued in [36], this length scale loses its 
supposed physical relevance. The former definition, based on the local rate of 
shear, is associated with the energy dissipation instead of with the enstrophy 
dissipation, thereby not suffering from this difficulty Nevertheless, given that 
there exists no dissipation length scale (for both energy and enstrophy) well 
separated from the forcing scale for a broad class of forcing [36] , including the 
one considered by the classical theory of turbulence [20,22,29], the terminology 
"dissipation length scale" may not refer to where in wavenumber space the 
dissipation of energy (or enstrophy) occurs. Rather, it refers, in a qualitative 
sense, to the length scale below which viscous effects dominate nonlinearity. 
We will see later in this section that this length scale is better quantified by 
the concept of determining nodes [13,14]. In fact, it can be identified with a 
"minimum length scale" of the flow that can be rigorously deduced. 

The concept of correlation length has a relatively long history in fluid dynamics 
(dating back to early in the 1900s), and in dynamical systems theory in general. 
There has been more than one definition of this length scale; two that are most 
relevant to the present discussion will be briefly reviewed. The first and simpler 
one is the two-point correlation length defined in terms of the correlation 
function 



C(x,x') EE ((0OM) - (cf>) x )((f>(x',t) - (Cf>) x )) t , 

where the angle brackets with subscript x denote the spatial average, the angle 
brackets with subscript t denote the time average, and (f>(x, i) is some local 
variable. In cases where C(x, x') ~ exp{ — \x — x'\/£} as \x — x'\ — > oo one can 
then define the two-point correlation length £. This length scale, if it exists, 
may be identified with the size of weakly-interacting subsystems. The second 
correlation length is called the dimension correlation length, the definition of 
which is due to Cross and Hohenberg [8]. According to the idea of extensive 
chaos, there exists a dimension density S being the ratio of the attractor fractal 
dimension (or equivalently Hausdorff dimension) D to the system volume V, 



Since 5 has the physical units of inverse volume, Cross and Hohenberg 
suggest defining a dimension correlation length £5: 
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6 = <rv- 



(12) 



where d is the dimension of the system's physical space. It is highly likely 
that the dimension correlation length so defined might not be a new length 
scale. The justification for this claim lies in the heuristic argument in favor 
of extensive chaos provided in Section 1. It is argued there that the attractor 
dimension D is of the order of (L/£) d , so that 5 = £~ d . Substituting this into 
(12) one obtains £5 = I. Thus, def. (12) might just be the familiar dissipation 
length scale. 

The present problem basically has four physical parameters, which can be 
independently controlled: the system linear length scale L, the forcing length 
scale £ s , the forcing root mean square amplitude (forcing density) |/| = ||/|| /L, 
and the molecular viscosity v. The first two parameters make obvious and nat- 
ural length scales. The last two parameters can be grouped, on dimensional 
grounds, to form a third length scale. Since the quantity iy 1 l\j\) x ^ has the 
physical units of length, it is tempting to designate this quantity as a dissi- 
pation length scale. It turns out, however, that by grouping u, |/|, and £ s one 
can form a more physically plausible length scale which will be identified with 
the dissipation length scale £: 



( z/ 2 \ 1/2 

Hwj • (13) 

This agrees with the definition of £ in [16] because, according to (9), \f\£ s /u 
is an estimate of (strictly, an upper bound for) the average rate of shear 
(||-u|| 2 ) 1//2 /L, thereby related to the energy dissipation (not enstrophy dissipa- 
tion). Thus, this definition is expected to differ from the one due to Kraichnan 
and Batchelor. To see the extent of this difference, we write out the explicit 
expressions for these length scales, denoting the one due to Kraichnan and 
Batchelor by £ KB and preserving the notation £ for the other, 



vL ^ 1/3 



TIT^T/? , (14) 



,x„0 1/2 / 



It follows that 



3 KB „ (\\<) 1/2 



(16) 
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The right-hand side of (16) is the enstrophy dissipation length scale, which ac- 
cording to the classical theory [3,20] would presumably be £rb- This would im- 
ply £rb ~ £■ However, in the usual case of a spectrally localized forcing around 
a forcing length scale £ s , such that the energy and enstrophy injection rates are 
related by (J n ufdx) = ^(/ n (-A)«/da;), one has (Nli> V 7(Nl2> 1/2 = U 
[36]. This result implies that £kb ~ (£ 2 £ s ) 1 ^ 3 , which makes the difference be- 
tween £kb and £ significant if £ <C £ s . 

Closely related to the notion of dissipation length scale are the concepts of 
determining modes, nodes, and finite-volume elements [11,13,14,18,19], from 
which a "minimum length scale" of the flow can be rigorously deduced. Here we 
restrict our discussion to the concept of determining nodes, which is introduced 
by Foias and Temam [13] and for which a sharp (sharpest to date) estimate of 
such a length scale is derived by Friz and Robinson [14]. A simplified version 
of Friz and Robinson's result can be stated as follows. Let {xi}^ =1 be a set of n 
lattice sites (nodes) equally spaced in the domain Q, so that the lattice spacing 
is L/n 1 / 2 . Suppose that u and v are two trajectories on the global attractor 
and that u and v agree at the lattice sites, i.e. u(xi) = v(xi), for i — 1, 2, • • • , n. 
If this condition is satisfied, where n is of the order of the attractor dimension 
D and the forcing is analytic, then u = v. This result means that a trajectory 
on the global attractor is completely determined by its values at D discrete 
nodes, hence the term "determining nodes". Obviously, the lattice spacing 
LjD 1 ! 2 coincides with the dimension correlation length scale £5 previously 
mentioned. 

In terms of the length scales described in the preceding paragraphs, the results 
of the last section can be written as 



D ^\w) ( 1+ln «r) • < 18 > 

where a constant of order unity has been dropped from each expression. A 
couple of remarks are in order. First, not only does expression (17) agree 
with the idea of extensive chaos, but it also takes the form suggested by the 
heuristic argument in the introductory section. Second, it is clear that the 
logarithmic correction term in expression (18) surely diverges as L — > 00; 
however, (18) may well be optimal for small domain size. In that case the sys- 
tem is slightly "super-extensive". When L is sufficiently large, the bound on 
D behaves extensively as (17) becomes optimal. In this regime, the attractor- 
dimension estimate grows linearly with the domain area, provided £ is fixed. 
Constancy of £ in a varying domain size requires fixed |/|, v, and £ s . A fixed 
£ s means that the non-dimensional characteristic wavenumber of the forcing 
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increases in exact proportion with L. Note that for the extensive regime (17) 
implies that £ < L/D 1 ^ 2 , so £ can be identified with the "minimum length 
scale" in the sense of Foias and Temam and of Friz and Robinson. For the 
super-extensive regime £kb plays the role of this length scale if the logarith- 
mic correction is ignored. The correspondence between £kb and the rigorous 
"minimum length scale" (ignoring the logarithmic correction) has previously 
been noted by Friz and Robinson [14]. The length scale that separates the two 
regimes may be called the extensivity length scale [34]. Letting L e denote this 
length, a comparison of (17) and (18) gives 



T \ 2 ( f3 \ 2/3 / f3 X 1 / 3 



t) \£ 2 £ s \ £H S 



It follows that 



Le ~ eXP l#p 4)1/3 ' (19) 

provided £ <C £ s . It can be confirmed that L e > £ s , so a non-trivial extensivity 
length scale does exist. 

Finally, the dimension correlation length as defined by Cross and Hohenberg 
[8] and advanced by Egolf [9,10] is just the newly defined dissipation length 
£, which has the physical significance as discussed above. In particular, £ rep- 
resents the spacing of determining nodes, and in the absence of a dissipation 
range well separated from the forcing region [36], the terminology associated 
with £ finds itself a better justification in this new sense. 



4 Conclusion 



We have analysed the two standard attractor-dimension estimates of the two- 
dimensional Navier-Stokes equations, derived by Constantin-Foias-Temam. 
It is shown that one of the estimates grows linearly with the domain area, for 
sufficiently large systems, when the kinematic viscosity and forcing density 
and its length scale are held fixed. This is consistent with the central idea of 
extensive chaos: namely the linear scaling of complexity, as characterized by 
the attractor dimension, with domain area. This qualifies (10) as the optimal 
estimate for large systems. For small systems, for which the domain length 
scale is below a threshold called the extensivity length scale, the estimate (11) 
becomes optimal. In this regime, the attractor dimension is slightly super- 
extensive, essentially by a logarithmic factor. 
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In the extensive regime, the upper bound on the attractor dimension is given 
by the ratio of the domain area to the square of the dissipation length scale, 
which is defined, on physical grounds, in terms of the average shear. This im- 
plies an equivalence between the dissipation length scale and the dimension 
correlation length scale. The latter has been suggested to be a new character- 
istic length scale for extensively chaotic systems, but it seems likely that it is 
just another way of determining the dissipation length scale. 
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